Asymmetric Fluid Criticality I: Scaling with Pressure Mixing 
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The thermodynamic behavior of a fluid near a vapor-liquid and, hence, asymmetric critical point is 
discussed within a general "complete" scaling theory incorporating pressure mixing in the nonlinear 
scaling flelds as well as corrections to scaling. This theory allows for a Yang- Yang anomaly in 
which /i^(r), the second temperature derivative of the chemical potential along the phase boundary, 
diverges like the specific heat when T Tc, it also generates a leading singular term, \t\^^ , in the 
coexistence curve diameter, where t = {T — Tc)/Tc. The behavior of various special loci, such as 
the critical isochore, the critical isotherm, the fc-inflection loci, on which x'*"' = x(P,r)/p*= (with 
X = p^k-QTKr) and Cy'' = Cv{p,T)/ p'' are maximal at fixed T, is carefully elucidated. These 
results are useful for analyzing simulations and experiments, since particular, nonuniversal values 
of k specify loci that approach the critical density most rapidly and reflect the pressure-mixing 
coefficient. Concrete illustrations are presented for the hard-core square- well fluid and for the 
restricted primitive model electrolyte. For comparison, a discussion of the classical (or Landau) 
theory is presented briefly and various interesting loci are determined explicitly and illustrated 
quantitatively for a van der Waals fluid. 



I. INTRODUCTION AND OVERVIEW 

In 1964 Yang and Yang [1] derived the thermodynamic 
relation 



Cy{T, V) = VT 



{dry 



{dry' 



(1.1) 



for a fluid at pressure p, chemical potential and tem- 
perature T, where V is the volume and N the number of 
particles, while Cy * is the constant-volume specific heat 
or, better, heat capacity. This has since been called the 
Yang- Yang relation [2]. When it is applied to the two- 
phase region beneath the critical temperature T^, one has 
p = Pa{T) and /i = /^o-(T) (where a denotes the phase 
boundary on which liquid and vapor may coexist), and 
the partial derivatives become total derivatives. Since the 
observations of Voronel' and coworkers in 1962-63 [3] it 
has been well established that the heat capacity Cy ' (T) 
diverges weakly along the phase boundary when the crit- 
ical point is approached. The divergence of Cy * then 
implies that one, the other, or both of the second deriva- 
tives pJ^(T) and fi'l{T) must diverge when T ^ T^— along 
the phase boundary. The lattice gas model and its stan- 
dard variants predict that jjl^ remains finite at while 
diverges like the specific heat [4]. However, Yang and 
Yang suggested that in real fluids both should diverge 
[1]: clearly this is a basic issue for the description and 
understanding criticality in fluids. 

Recently, Fisher and coworkers [2] carefully analyzed 
experimental two-phase heat-capacity data for propane 
(CaHg) and CO2 and showed that the evidence rather 
strongly indicates that /i^ does indeed diverge like the 
specific heat when T — > T^—. They dubbed this phe- 
nomenon a Yang- Yang anomaly [2] . Even though im- 
purities in the propane system have definite effects on 



the heat-capacity data [5] , the existence of a Yang- Yang 
anomaly cannot be ruled out and , in our assessment, 
remains the most plausible scenario. In fact, Orkoulas et 
al. [6] have performed grand canonical Monte Carlo sim- 
ulations for the hard-core square- well fiuid and concluded 
that this model system probably exhibits a negative but 
small Yang- Yang anomaly, i.e., a specific-heat-like diver- 
gence in the chemical potential derivative, d^^^/dT'^, 
of magnitude significantly less than the divergence of 
v{d'^Pcr/dT^), where v — V/N . How can one then accom- 
modate such a Yang- Yang anomaly in scaling theory? 

The concept of asymptotic scaling has proved a pow- 
erful tool for gaining insight into critical phenomena in a 
variety of systems including fluids [7-12]. Furthermore, a 
scaling equation of state has been rather well conflrmed 
experimentally for many fluids [13,14]. The currently ac- 
cepted asymptotic scaling description of fluid criticality 
[15] requires two scaling fields, namely a thermal field, t, 
and an ordering field, h, that, in leading order, are both 
linear combinations of t (x T — and h = fi — fic^ 
the deviations of the temperature and chemical potential 
from their critical values. This description has also been 
extended to describe fluctuations in finite systems and 
applied to estimating the critical points of model fiuids 
[16]. However, within this scaling description, the chemi- 
cal potential is always analytic along the phase boundary 
and through the critical point, thereby having a finite 
value of the second temperature derivative, fi'^. Hence, 
in order to account for a Yang- Yang anomaly, the current 
scaling description must be modified. Indeed, Fisher and 
Orkoulas [2] argued that the pressure deviation, p — pc, 
must mix into the scaling fields, especially into the or- 
dering field h. 

In order to formulate a scaling theory in which the 
pressure is mixed into the scaling fields, let us consider a 
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full therm,odynamic description of a one-component fluid 
as provided by a functional relation between the three 
thermodynamic fields, pressure, p, chemical potential, fi, 
and temperature, T, say 



T(p,M,T) = 0. 



(1.2) 



Extending the approach sketched in Ref. [2] to include the 

full spectrum of correction exponents (see also [17,18]) 
one may formulate a rather general scaling hypothesis 
by asserting that near a typical critical point, (pc, Mc, 7c), 
the thermodynamics can be described, at least asymptot- 
ically, by 

4'(A2-"p, At, A^/i; X-^^U4, X^^'U5, ■ ■ ■) = 0, (1.3) 

where A is a free, positive scaling parameter. In this ex- 
pression p{p, fi, T), i{p, fi, T) and h{p, /i, T) are the three 
relevant nonlinear scaling fields, while u^lp, ii,T) and 
U5{p,iJ,,T) are the leading even and odd irrelevant scal- 
ing fields, respectively. Including the quadratic nonlinear 
terms, we may write the basic nonlinear fields as 



p = p~ kot- k) fi 

- rat^ - qa fi^ - votfi - mop^ - nopt - nspfi 



+ 03{t,ii,p), 
i=t-lifi- jip 



(1.4) 



rit^ — qijl^ — vitfi — m\p^ — nipt — nipfi 



+ Oi{t,jl,p), 
h = p-kit- j2p 



(1.5) 



r2t^ - 92^^ - '>^2til — m2p^ — n2pt - n^pfi 



(1.6) 



where we have introduced the dimensionless critical de- 
viations [2] 
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T-Tc 



IJ.= 



P = 



P-Pc 
Pc^bT'c 



(1.7) 



in which pc is the critical (number) density and ks is 
Boltzmann's constant, while Om{x,y,z) denotes a for- 
mal expansion in powers x^y^z^ with 7 + fc + Z > m. In 
accepting these expansions we are neglecting any scaling- 
exponcnt "resonances" that might complicate the expres- 
sions by introducing logarithmic factors, etc. [19-22]. For 
(d=3)-dimensional systems, particularly real asymmetric 
fluids — which are of especial interest — this should be 
quite satisfactory. The irrelevant scaling fields, W4, U5, 
• • •, will, in general, have similar expansions in powers of 
the variables, fi and p. 

As usual, a in (1.3) is the universal critical exponent 
of the specific heat and A is the gap exponent that is 
related to a and the exponents [3 (for the spontaneous 
magnetization), 7 (for the susceptibility) and 5 (for the 
critical isotherm) by the scaling relations 



2-a-/3 = /3 + 7 = /3J, 



(1.8) 



while 04^ = and ^5 arc the positive leading corrcction-to- 
scaling exponents. In the case of the ((i=3)-dimensional 
Ising universality class, which is believed to character- 
ize fluid criticality rather generally, we may accept a ~ 
0.109, P ~ 0.326, 7 ~ 1.239, A ~ 1.565, ~ 0.52 and 
05 ^ 1.32 [17]. Finally, substituting A = into (1.3) 
yields the "thermal scaling" form 

$±(p/|iT-«, h/\i\^; u4if, W5|iT=, • • •) = 0, (1.9) 

where $^(a;, y, ■ ■ ■) = '^{x, ±1, y; • • •) in which ± corre- 
sponds to 0. 

A further observation is worth mentioning at this 
point. Specifically, if the critical point is drawn out into 
a lambda line without changing the universality class, 
by, e.g., the imposition of a magnetic field, H, or, the 
addition of a second molecular component with chemical 
potential, say p}^, etc., the only change needed in the for- 
mulation is the inclusion of a further, nonordering field 
g [ (X H 01 zu = ex.p{—fj}^/kBT), etc.] in the expansions 
of the three nonlinear scaling fields; likewise for any fur- 
ther fields that leave the universality class unaltered: see, 
e.g., [17,18]. 

Now the particle number N, volume V, and the en- 
tropy S, are related by the Gibbs-Duhem relation 



Vdp - SdT - Ndp. = 0. 



(1.10) 



Hence the number density and the entropy density are 
given by 



P = 



N 
V 



dp^ 



S = 



V 



dp 
df 



(1.11) 



Other thermodynamic quantities, such as the specific 
heat, compressibility, etc., follow in the standard way. 

The crucial point about the scaling formulation pre- 
sented here is that the three standard thermodynamic 
fields, p, /i, and T enter in a fully symmetrical way with 
no a priori assumptions as to which pair combination 
or thermodynamic potential, say p{j2,T) or ii{p,T), is 
"more basic." By contrast, previous formulations have 
typically concluded that it was most appropriate to re- 
gard p(/x,T) as the "basic" quantity [15,17]. This is 
what one is led to by studying the standard lattice gases 
and what arises most naturally from field-theoretic and 
rcnormalization group approaches [20-22]. However, one 
might note that the original Widom formulation [7] was 
based on integrating the equation of state and led, fairly 
naturally, to a scaling hypothesis for the Helmholtz free 
energy density /(/?, T); on the other hand, the exactly 
soluble cluster-interaction fiuids of Fisher and Felderhof 
[23] showed that fi{p,T) was the appropriate thermo- 
dynamic potential for describing critical-point scaling in 
these models. 

It is easily seen, however, that the previous p(/i, T) for- 
mulation [15,17] is recaptured here simply by suppress- 
ing the p dependence in the nonlinear scaling fields h 
and i: i.e., in linear order, by setting ji = j2 = 
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in (1.5) and (1.6), and in quadratic order, by setting 
mo ~ mi — ■ ■ ■ = nr^ = 0. In that case the expansion of 
p in (1.4) merely serves to represent the usual "smooth 
background" po{n,T) that is always required. 

The original discussion of scaling in fluids [7,8] had 
to incorporate in the definition of the ordering field h 
the mixing coefficient ki which, indeed, is then directly 
proportional to {diJa/dT)c where ndT) represents the 
phase boundary below T^,; but the coefficients ji, j2 and 
h did not appear. Studies of various models later re- 
vealed that the coefficient Zi, which mixes the chemical 
potential into the thermal scaling fields t, should also 
be included: doing so yielded an unexpected singular- 
ity in the coexistence curve diameter, p(T), proportional 
to (See, e.g. Ref. [15].) Nevertheless, it was still 

argued (in Ref. [15b]) that one could suppress the linear 
pressure mixing coefficients ji and j2 despite the exis- 
tence of soluble models in which scaling, in fact, required 
them [23]. 

It is only the recent reconsideration of the possibility 
of a Yang- Yang anomaly in real fluids and nonsymmetric 
models [2,5] that has led to the realization that pressure 
mixing, i.e., nonvanishing cocfhcicnts ji and j2 in (1.5) 
and (1.6), should be reconsidered. Indeed, as shown in 
Ref. [2], if a nonvanishing Yang- Yang anomaly arises, so 
that ((i^/i(j/(iT^) diverges when T — > Tc — , then, within a 
scaling formulation, one must have j2 7^ 0. In fact, that 
also leads, as we demonstrate below, to a singular term 
varying as \t\^^ in the coexistence curve diameter, p(T), 
which, in realistic situations, will actually dominate the 
previously discovered Ifl""^"" term [2,15]. Conversely, the 
absence of pressure mixing in the linear ordering field, i.e. 
j2 = 0, implies the absence of a Yang- Yang anomaly as is 
the case for the standard lattice gases. Granted that i2 
may not vanish, it is clearly appropriate to allow ji ^ 0, 
i.e., to consider pressure mixing also in the thermal field, 
i. 

In light of this background the aim of the analysis pre- 
sented in this paper is to thoroughly explore the implica- 
tions of the "complete scaling hypothesis" embodied in 
the relations (1.3)-(1.9). Hence in Sec. II the "complete 
scaling" theory is formulated with the incorporation of 
pressure mixing as well as corrections to scaling. For 
this purpose, we mostly follow the treatment given in 
Ref. [17]. The properties of the scaling functions in the 
two limiting cases, h/\t\^ or 00, are set out as 
well as the consequences of the thermodynamic convex- 
ity or the 2nd Law of Thermodynamics [24] . For use in 
the subsequent calculations, generalized scaling densities 
and susceptibilities are defined and the relations between 
them and the actual physical quantities are determined. 

Once the formulation of the scaling theory is com- 
pleted, we examine, in Sec. Ill, the standard singu- 
lar thermodynamic properties including the coexistence 
curve and the specific heat. We explicitly determine the 
nature of the phase boundaries, Pa-{T) and fia{T), of the 
coexistence curve diameter, p = ^(puq + Pvap) and of 



the half-jump. Ap = ^(puq — Pvap), as well as of en- 
tropies at the coexistence curve. As mentioned above, the 
pressure-mixing coefficient j2 generates a singular 
term in the phase boundary Pa{T) and thereby, the di- 
vergence of the {d^ Her / dT^) at the critical point: see Sec. 
III. A; in addition, a singular \t\'^^ term in the coexis- 
tence curve diameter appears and dominates the \t\^~°' 
term known previously [15]: see Sec. III.B. In Sec. III.C 
the linear mixing coefficients ji, j2, • • •, h are related to 
various physical amplitudes, etc., which can, in princi- 
ple, be measured via experiments or simulations. The 
Yang- Yang anomaly, that was the main motivation for 
this study, is discussed in Sec. III.D; it is shown that the 
strength of the anomaly, namely 7^^, which measures the 
contribution of the chemical potential to the heat capac- 
ity relative to that of the pressure in an asymptotic limit 
— see (1.1) — is related solely to the pressure-mixing 
coefficient j2- 

In Sec. IV we study a number of special loci, partic- 
ularly in the density-temperature (p, T) plane, that in- 
tersect the critical point. These loci are of interest in 
their own right but they prove to be especially useful in 
attempting to estimate precise values of the critical den- 
sity, pc, both in the case of real asymmetric fluids and in 
the simulation of various model systems such as the hard- 
core square-well fluid [6] and, more challengingly, primi- 
tive model electrolyte systems [25]. More concretely, we 
study the locus ^ = pc, dubbed the critical isokyme, the 
critical isobar {p = p^), the critical isotherm (T — Tc) 
and the critical isochore {p = pc) in the {p,T), {p,T), 
{p,T), etc., planes. The singular exponents characteriz- 
ing these loci are obtained and the corresponding ampli- 
tudes are expressed in terms of the mixing coefficients 
and the expansion coefficients of the scaling functions 
and, thereby, related to one another. One discovers, how- 
ever, that these characterizations of the various loci may 
not be helpful in estimating the critical point in practical 
applications to experiments or simulations, since various 
critical parameter values must be known a priori: but, 
except for exactly soluble models, the requisite values 
are not normally available. 

Hence, for direct applications to the analysis of sim- 
ulation data, we also study in Sec. IV. E, the "fc- 
inflection loci" introduced by Orkoulas et al. [6] . The "fc- 
susccptibility loci" are defined by the points of isothermal 
maxima of the modified susceptibility x*-*^-* = xl above 
Tc in the {p,T) plane where x = P^k-QTKT is the stan- 
dard isothermal susceptibility while Kt is the compress- 
ibility. We find that these loci have leading \t\'^^^ terms 
followed by more slowly varying and t contribu- 

tions: however, the leading amplitude vanishes when k 
takes some 'optimal value' fcopt- Thus, when k = kopt, 
the corresponding A:-locus "points" most directly to the 
critical point. Furthermore, we find that /copt is directly 
related to the strength of the Yang- Yang anomaly via 
^opt = 37?.^. With the aid of simulations we illustrate 
these loci quantitatively for the hard-core square-well 
fluid and the restricted primitive model electrolyte, and 
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compare the results with those for a van der Waals fluid: 
see Figs. 1-3, below. 

One may, similarly [26], define "/c-heat-capacity loci" 
via the points of isothermal maxima in the (p, T) plane 

of the modified specific heat Cy^ = Cy j . These loci 
have not been used in past simulations but may be useful 
in the future. We find, as for the susceptibility loci, that 
a leading \t^^ term appears in the near-critical expan- 
sions, followed by and t terms. Again, the leading 
amplitude vanishes when fc is equal to a special value that 
is once more related directly to the Yang- Yang anomaly: 
see Sec. IV.F. 

In order to provide semiquantitative guides to the be- 
havior of these loci in real systems, we study them in Sec. 
V, using classical theory for the gas-liqnid phase transi- 
tion on the basis of a Landau order-parameter expansion 
of the free energy. We obtain the vapor pressure, Pa{T), 
and the saturation chemical potential, i.i^{T), up to sec- 
ond order in t, while in the current literature they appear 
only up to first order [13]. The coexistence curve diame- 
ter, p{T), is also found to order t. These curves can be an- 
alytically continued from T < Tc to the one-phase region 
above Tc so providing special loci that cannot, in general, 
be defined in nonclassical cases. Using the van der Waals 
equation, these analytic extensions are illustrated quan- 
titatively: sec Figs. 5-7. In addition to these critical loci, 
the /c-susceptibility-loci are obtained and compared with 
those of the hard-core square-well [6] and the restricted 
primitive model [25]. In Sec. V.D the Yang- Yang relation 
(1.1) is discussed in more detail within classical theory 
and extended to general loci in the {p, T) plane. Finally, 
a Yang- Yang-type relation along the critical isotherm is 
discussed briefly. 

Sec. VI summarizes the paper and provides a key to 
the main results. Some further details are provided in the 
Appendix and others are available in the first author's 
thesis [26]. 

II. SCALING FORMULATION AND 
THERMODYNAMIC FUNCTIONS 

To explore the scaling description (1.9), it is appropri- 
ate to focus on the relevant scaling variable (or combina- 
tion) 

y{p,p,t) = Uh/\i\^, (2.1) 

where, without loss of generality, U may be taken as a 
positive constant. The basic reason for this choice, as 
against p/\t\'^~°', is that A is, in general, less than (2 — a) 
[since, see (1.8), (3 > 0] so that y diverges less rapidly 
when t — > 0. Beyond y we need to account for the full 
set of irrelevant scaled variables, namely, 

yfe(p,/i,^) = ^7fc(p,A,^)|^T^ 

0k+i>Ok>O, fc = 4,5,---, (2.2) 



where Uk oc Uk and we will assume, when needed, that 
the associated irrelevant amplitudes Uk are noncritical, 
meaning that Uk{p,fi,t) can be expanded in a formal 
power series of p, fi and t (which may not necessarily 
converge). Since a fluid has, in general, no obvious sym- 
metry, we do not impose any restriction on the Uk'- but 
see also the discussion in Ref. [17] Sec. II. 

Now the thermodynamic potential for the fluid (in this 
case the pressure — as a consequence of selecting y as 
the primary scaled variable) can be written by formally 
solving (1.9) as 

P = Q|i?~"T4^±(y,2/4,y5,---), (2.3) 

where Q is a positive amplitude while W± is a scaling 
function that, in the case of bulk fluids, embodies the 
properties of the {d = 3)-dimensional Ising universality 
class. In this expression the subscript ± refers to f ^ 0. 
On choosing appropriate values for Q and U , the scal- 
ing function, W±{y,---) becomes universal. Note that 
the usual analytic background part of the potential, say 
Po{t:P'), is included in the nonlinear scaling field p: see 
(1.4) which may, clearly, be solved iteratively for p to 
yield an expansion for po{t, fi). 



A. Scaling functions 

The scaling function W^±(2/, j/4, ys, • • ■) should be both 
universal and invariant under change of sign of the odd 
arguments y, y^, yr, Since the irrelevant scaling 

variables, 7/4, 1/5, • • •, vanish as the critical point is ap- 
proached, we can also expand the scaling function W± as 
[17] 

W±{y,y4,y5, ■ ■ •) 

= (y) + yiWi^^ {y) + 1/5 W^i'^ iv) + ■ ■ ■ 

+ ylwi''^\y) + y4y5Wi^''\y) + --- 
= ^M/f(2/)2/W, (2.4) 

where for brevity we have used the multi-index, k, de- 
fined via 

K = 0, (4), (5), • • • , (4, 4), (4, 5), • • • , (4, 4, 4), • • • , (2.5) 
as a label and as an exponent via 

yO^l, yli^j>-M=y.y....y^_ (2.6) 

We consider k = (i, j, • • • , n) to be odd or even according 
to whether the sum i + j + ■ ■ ■ + n is odd or even. The 
symmetry of W±{y, 2/4, 2/5, • • •) then requires [17] 

W^i-y) = i-rw^iy). (2.7) 

Recognizing this symmetry, one may write expansions 
for W^{y) for small y. For f > we have 
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= + y^W^^ + y^H/+4 + • • • , for /« even, 

= yW^+i + + + • • • , for k odd. (2.8) 

By choosing appropriate values for the nonuniversal met- 
ric ampUtudes Q, U ^ etc., the expansion coefficients can 
be normaUzed, in general, such that [17,24] 

W°2 = W'+o = 1 i.>^ even) or W^-^ = 1 (k odd). 

(2.9) 

For i < 0, the existence of the first-order transition 
leads to \y\ factors in the expansions so that one has [17] 

wiiy) = [W^o + Iviw^i + + ly|'w^-3 + • • Wiv), 

(2.10) 

where the special signum function is defined by 



o-«(y) = 1 

= sgn(?y) 



for K, even, 
for K odd. 



(2.11) 



Thermodynamic convexity (which embodies the Second 
Law) then requires that W^-^ and W% must both be 
positive: sec Rcf. [24]. 

For large arguments, \y\ oo, the scaling functions 
W!^{y) and W!l{y) must satisfy stringent matching con- 
ditions to ensure the analyticity of the potential through 
the plane t = for all ft, ^ 0. By these conditions, the 
functions W^{y) may be written as 



i+>;(±)'«;riy|-'/^ 



1=1 



cTniv), (2.12) 



where the multiexponent 9[k\ is defined by 



^[0]=0, 0[{i,j, 



, n)] =ei+0i + 



(2.13) 



with 



^, J, 



■ , n > 4. By virtue of the normalizations 



(2.9) the numerical amplitudes W!:j, W!^^, W^, and 
should all be universal. Beyond that, thermodynamic 
convexity dictates that and ^2 must be positive 

while {wi)^/w2 must be bounded above [24]. The sign of 
Wi is not determined by convexity alone but must, in gen- 
eral, be negative: see Ref. [24]. This plays an important 
role in determining allowable phase diagrams in a density 
(or composition) space at critical endpoints [17,24]. 



B. Generalized density and entropy 

To extract explicit results from (2.3) revealing the sin- 
gularities of the various thermodynamic derivatives, crit- 
ical loci, etc., we first recall (1.11) and define a dimen- 
sionless reduced density, p, and entropy, s, by 



P = 



dp 
dfi 



Pck^ 



(2.14) 



Similarly, the generalized number density, p, and entropy 
density, s, will be defined by 



P = 



dh/ t 



(2.15) 



These quantities then obey the usual simple scaling rules, 
namely, p ^ jij'^ and s ^ when f ^ 0. 

To find concrete expressions for p and s, one may con- 
sider the differential relation. 



dp = pdh + Mi, 



(2.16) 



where it is appropriate to recall that all the scaling vari- 
ables yk in (2.3) are functions (only) of p, p, and t and, 
hence, of p, p and t. Using (1.4)-(1.6), this differential 
relation may be written as 

(1 - 2mop - nop - n^t + j2P + jis H )dp 

= {h + nop + 2qop. + vot + p - hs -\ )dp. 

+ {konsp + vop + 2rot + s-kip+-- ■)dt, (2.17) 

where in the brackets we have retained only terms up to 
linear order. Using (2.14), we then obtain the relations 



(2.18) 
(2.19) 



, _ Zo -I- nop + 2qop + Vot + p - hs + ■ ■ ■ 

^ 1 - 2mop - nop - nst + j2P + jis -\ 

ko + n3p + vop + 2rot + s - kip 

s = 

1 - 2mop - nop - n^t + j2p + jis -\ 

If the nonlinear mixing coefficients m,j, nj, • • •, Vj in 
(1.4)-(1.6) are ignored, these expressions may be approx- 
imated by 



lo + p- hs 
1 + j2P + jis' 



ko + S - kip 
1 +i2/0 + iis' 



(2.20) 



Note that p and s are nonlinear functions of p and s 
owing to the pressure mixing coefficients, ji and j2- In 
the absence of pressure mixing {i.e., ji = j2 = 0), we 
recover the Bruce- Wilding linear relations [16]. Wo will 
see later that the approximation (2.20) is adequate to 
derive the most important singularities in leading order. 



C. Generalized susceptibilities 

The second derivatives of the potential determine the 
usual response functions, e.g., susceptibilities, heat ca- 
pacities, etc. We define the basic reduced susceptibilites 

XNN = {d'^p/dp?)t, xuu = {d'^p/dt'^)jx, 

XNU = {dH/dpdt) = {d^p/dtdfi). (2.21) 

These are related to the number and energy fluctuations 
most directly accessible in grand canonical simulations 
[6]. The isothermal susceptibility x is defined by 
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(2.22) 



III. THERMODYNAMIC PROPERTIES 



and is related to the isothermal compressibility, Kt = 
p~^{dp/dp)T, by X = P^Kt and to the reduced suscep- 
tibility Xnn by 



X = iPc/kBTc)XNN- 



(2.23) 



Experimentally, the most important heat capacity for 
fluids is the constant- volume heat capacity with a density 
defined by 



T fdS 



(2.24) 



where S is the entropy density defined in (1.11). It is 
convenient to introduce the dimensionless reduced spe- 
cific heat, namely. 



Cv = {ds/dt)p, 

for which one has 

pCv/T = {kBPc/n)Cv. 



(2.25) 



(2.26) 



The reduced specific heat, Cy, is then related to the re- 
duced susceptibilities in (2.21) via [6] 

Cv = Xuu - X%u/XNN. (2.27) 
Finally, we define generalized (scaling) susceptibilities 



via 



fd-'pX fd''p\ 

Xhh = -c- , Xtt = T37 , = ' ' 



'2p 

\dhdi) 

(2.28) 



From (2.3) one finds Xhh ^ \i\ ^' , Xht ^ \ and 
Xtt 1^1 ~" when t — !■ 0. When the nonlinear mixing 
terms in the scaling fields are ignored, as above, XAfiv 
can be expressed in terms of the generalized densities 
and susceptibilities as 

XiVJV « [(ei + e2s)'^Xhh + (es + e2pYxu - 2(ei + e2s) 
x(e3 + e2P)Xht] /(I + .72/5 + ,7lS)^ (2.29) 

where the derived mixing coefficients are 

ei = 1 - hh, 62 = ji + nh, es = h+ hh- (2.30) 

The detailed derivation of (2.29) is presented in Appendix 
F of Ref. [26] henceforth to be denoted K; there it is also 
shown how (2.29) [K(3.41)] captures the leading singular 
correction terms. 



A. Phcise boundaries 

The phase boundaries, say Pa{T) and Pa{T) [or Pa{t) 
and Pa{t)], on which two phases may coexist, can be de- 
termined by equating the pressure and chemical potential 
on the vapor and liquid sides below Tc {i.e., for i < 0). 
On using the small y expansion of p for t < with the 
aid of (2.10), we obtain (including higher order terms) 

P± = Q\i\'-'iw% ± W'l.y + w°2y^ + ■■■ 

± C/5c|if ^ {W['^ ± W[% + W['^y' + ...) + ...], (3.1) 

where ± now refers to h ^ 0, while U^c and U^c arc the 
critical values of the irrelevant scaling amplitudes, U4 and 
U5, respectively. Equating p+ and p- then yields 



(3.2) 



where one sees that various even and odd terms cancel. 
Solving for y iteratively, one obtains 



(5) 



+ • 



(3.3) 



Note that in contrast to the symmetric case (e.g., the fer- 
romagnetic Ising model) where U5 = 0, the scaling field 
does not vanish identically along the phase boundary. 
Using the definition (2.1), we find the phase boundary is 
given by 



K{t) = Ei\t. 



z |A+e5 



7 lA+e+es 



where the coefficients are 



a =- 



Eo = 



J7(Vl^_i 



\2 



(3.4) 



(3.5) 



Now substituting p, — Pa{T) and p = PaiT) in (1.5) 
and (1.6) and using (3.4), we obtain an expansion [see 
K(3.65)] of 



Ma(T) = [/i<.(T)-/ie]ABTe 

in powers of t, of pa itself, and of 

P<y{T) = [Pa{T) - Pc]l PcknTc 



(3.6) 



(3.7) 



One can solve this equation iteratively for p^ as a func- 
tion of t and Pa to obtain 
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iJ-TiT) = j2pa + kit+{m2+ 3lq2 + j2n2)pl 
+ {r2 + fci92 + V2ki)t^ 
+ (ns + 2j2kiq2 + kin2 + V2j2)Pat 



+ Ei\{l-kih)t-{3i+j2ll)pa\ 
+ E2 I (1 - kih)t - {ji + j2h)Pa I 



Substituting this result into (1.4) yields 

P<7 = (1 - 32k)Pa - (ko + kilo)t 
- EiPl - Eit" - EaP^t + ■ ■ 



(3i 



(3.9) 



where the coefRcients Ej depend on the quadratic 
nonlinear-field coefficients ro, mo, •••in (1.4) to (1.6): 
see K(3.68)-(3.70). 

Similarly, from (1.5) and (3.8), one finds 



ta = {l- kih)t - {ji +j2h)Pa + 



(3.10) 



Substituting this into (3.1) with the aid of (2.1) and 
(3.4), combining it with (3.9) and solving iteratively for 
pc finally yields the result 



p^{T) = p^^it + pa,2t'^ + Ap\t\^-" [1 + ap\tf + ■ 



+ bp\tf 
+ ■■■, 



+ o.p\t\ 



65-13+6 



+ 



where the leading coefficients are given by 

E2.Pa,l, 



kg + klip _ F -2 , 



1 - .72^0 ' 



, ,2-a 



1 - J2^0 

y^-j2lo)£h^,e,. 



a„ = 



■/3 



QW% 

in which we have introduced mixing factor 

T = 1 - klh -Pa,l{jl +j2h), 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



while El is given in (3.5) and and bp can be derived 
from K(3.74). 

Note that the leading even correction-to-scaling term 
enters the phase boundary Pa{T) with an exponent (2 — 
a) + 6, while the odd correction-to-scahng term has an 
exponent (2 — a) — /3-|-^?5. In addition, it is not hard to see 
that the subsequent term 6p|tp^°^ (and analogous terms 
below) must be preceded by lower order terms such as 
dp\tr, a'X+\ ep\t\^<^, 4'\t\o+\ etc. 

Now let us substitute the result (3.11) back into (3.8). 
With a little further effort, one obtains the desired re- 
sult for the chemical potential on the phase boundary, 
namely, 



ii„{T) = fi^^it + /t<,,2i' + A^|t|2-«[l + a^|t|^ + • • • 

+ hX'~^ + a'X'-^^' + ■■■ + b'X'" + • • •] 
+ •••, (3.15) 

where the leading coefficients satisfy 

fl-a,! = kl +j2pa,l, = j2^p, = Up, (3.16) 

6^ = ^i|t|^+^V^^, a'^ = b^,cp, b'^ = b'p, (3.17) 

while jla,2 is given in K(3.78). 

As is to be expected, the spectrum of singular terms 
that appears in the expansion for /icr(r) is the same as 
that for pu{T). Note, however, that the leading singular 
amplitude, for HaiT) is proportional to the pressure- 
mixing coefficient j2 • Thus, in contrast to the traditional 
scaling treatment, pressure mixing in the scaling fields 
implies that the second derivative of Ha{T) diverges at 
the critical point with the same exponent a as the specific 
heat. However, even in the absence of pressure mixing 
per se the phase boundary ^a{T) is not analytic: rather 
its third derivative diverges (in contrast again to lattice- 
gas models) owing to the odd correction-to-scaling term, 
since in the case of fluids one has 2 < 2 — a — 13 + 0^ = A-)- 
6*5 < 3. [Note that the amplitude, of |t|2-a-/5+e5 

in (3.15) does not vanish when j2 = 0.] The conclusion 
that fi'liT) exhibits a cusp-like behavior near the critical 
point was originally advanced by Lcy-Koo and Green [27] 
and enters into the analysis of the effects of impurities 
on the detection of the Yang- Yang anomaly [5]: see the 
discussion below in Sec. III.D. 



B. Densities and entropies at coexistence 

The generalized densities along the phase boundary, po- 
and Scr, can be obtained from (2.3) and (2.15) by using 
the results (3.11) and (3.15) for p^ and jla- After some 
algebra, we find 

p^it) « ±QU\Ttf [wli + U4.w[^l\Ttf 

{w^^l-{wl2w['^/w'_i)} 



±u. 



5c • 



Tt 



(3.18) 



Sa{t) « -g|rt|i^" [(2 - a)W% + (2 - a + 9)Ui^W^^^\Tt\ 



±{A + 95)U5cW[''^\Ttf' + ■ 



(3.19) 



where ± refers to /i ^ 0. Using (2.18) and (2.19) and 

expanding in powers of t, we finally obtain the number 
density and entropy density on the two sides of the co- 
existence curve as 



|l-a 



+ Ait + 



+ A5\tf+^' + ■■■ 

± B\tf [1 + be\tf + b20\tf^ + ■■■]}, (3^20) 
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5±(T) = pcfcB{fco + 52/3 + + 5ii + • • • 

± B,\tf [1 + + d2M^^ + •••]}' (3-21) 

where terms varying as \t\^^^ and have not 

been displayed and higher order terms such as ±|i|''+'^, 
|^|/3+2-a^ etc. will also be present in general. Implicit in 
these results is the conclusion 

Iq= Pc = l and ko = Sc= S^/ Pck^- (3.22) 
The leading coefficients for the density are then 

B = {l-j2)QUWWTf, be = U^cW^^M'/Wl,, 

(3.23) 

b20 = j|sV(l - hf, A20 = - i2BV(l - h), 

(3.24) 

= (2 - a){h+n)QW%\T\'-'-, (3.25) 
Ai=vo+mz + {2qo + no)fia,i + (no + 2mo)pcr,i, 

(3.26) 

while the coefficients for the entropy satisfy 

5,/S = -(fci+j2fco)/(l-j2), de^bg, (3.27) 
d2i3 = jlBl/{k^ +j2kof, S20 = -32Bl/{ki +j2ko), 

(3.28) 

= (2 - a){hko - 1)QW%\t\'-'', 55 = B.A^/B, 

(3.29) 

and A5 and S5 follow from K(3.83, 3.84). 

From (3.20) the coexistence curve diameter, p{T) = 
5[Aiq(r) + /9vap(r)], and the width, 2Ap(T) = puq(r) - 
Pvap(r) of the coexistence curve, follows immediately. As 
anticipated, wc see that the diameter contains a ^"^^ 
term that is proportional to the pressure-mixing coeffi- 
cient j2'- and, since 2/3 < 1 — a for typical fluids, this ac- 
tually dominates the previously anticipated term 
[15]. Likewise, one may read off the entropy diameter, 
S{T), and entropy jump, AS(T), from (3.21). Again one 
observes that the dominating jfp'^ term in S{T) is pro- 
portional to j2- 

The total entropy in the two-phase region is given by 



ST\T) = V^-N^ 



(3.30) 



On the critical isochore, p = Pc, this yields the entropy 
density 

S^{T;p,)=Sc + PckB{l-j2)Ap\t\'-'' + Oit) (3.31) 

from which, as was to be anticipated, the terms have 
cancelled. [Recall that Ap is defined in (3.13).] 

In Fig. 1 grand canonical simulation data for the co- 
existence curve of the hard-core square-well fluid are 
presented. Adopting the critical point estimates, Tc ^ 



1.2179 and pc ~ 0.3067 [6], and the Ising values for 
the critical exponents yields the estimates B = I.20264, 
A2/3 = -0.00073, -4i_„ = 0.1897, Ai = -O.O69I4, 
bg = —0.2576, and 62/3 = —0.0852 for the amplitudes in 
(3.20) [6]. This fit is shown in Fig. 1 as a solid line that 
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FIG. 1. Coexistence data (open circles, with closed cir- 
cles for the diameter) for the hard-core square-well fluid ob- 
tained by Orkoulas et al. [6]. The solid lines connecting 
the data points represent an Ising-type fit: see Eq. (3.20). 
The vertical and horizontal dashed-lines locate the critical 
isochore, p* = pc — 0.3067, and the critical isotherm, 
T* = Tc ~ 1.2179 [6]. The curves above criticality and in 
the inset depict estimates for the fc-susceptibility loci for the 
values of k indicated. 

connects the coexisting density estimates (circles) to the 
critical point. Similarly, coexistence simulation data (in 
a C = 5 fine discretization level) for the restricted primi- 
tive model electrolyte are presented in Fig. 2. The solid 
line — the Ising fit to the coexistence data — is drawn 
by adopting the critical point values, Tc ~ 0.05069, 
Pc ~ 0.079 [25] and the amplitude estimates, B = 0.274, 
= 0.0165, = 0.919, Al = 0.586, bg = 1.464, 

and 62/j = 2.254. However, one must note that the spe- 
cific numbers attached to these particular amplitude es- 
timates cannot have a very significant meaning, unless 
the higher order corrections are considered more carefully 
than is practicable with the data currently available. 



C. Mixing coefficients 

It is clearly of interest to express the linear mixing co- 
efficients entering the scaling fields (1.4)-(1.6) in terms of 
various thermodynamic quantities which might, at least 
in principle, be measured in experiments or simulations. 
As regards the field p, we have already noted in (3.22) 
the simple expressions for fco and ^o- 

The pressure- mixing coefficient, j2, can be obtained 
via (3.16) from the singular amplitudes A^ in p,a{T) and 
Ap in Pair) simply as 
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FIG. 2. Plots of coexistence curves, fe-susceptibility loci, 
etc., as in Fig. 1, based on simulations for the restricted prim- 
itive model by Luijten et al. [25]. The critical parameters 
adopted are ~ 0.05069 and — 0.079 corresponding to a 
^ = 5 fine-discretization level [25]. 



j2 = Af,/Ap. 



(3.32) 



Then, from the observable Hmiting derivatives fi'^^ = 
{d^cr/dT)c and p'^^ = {dpa/dT)c, which correspond to 
ficr,i and pa.i, wc may obtain 



(3.33) 



The remaining two linear mixing coefficients, ji and li, 
can be obtained by using the ampUtudes, Ai^a and 5*1-0; 
in (3.25) and (3.29) which describe the singularity 
in the density and entropy diameters. Taking a ratio 
yields 



h+ji Ai_ 



jiko - 1 5i_ 



(3.34) 



where we have used = 1, while fco is given in (3.22). On 
the other hand, from the ratio of Ai 
via (3.13) and (3.25) the distinct relation 



to Ap we obtain 



- Ji 



Al. 



(2-a)(l-j2)Ap' 



where r is given in (3.14). Since r is linear in the mixing 
coefficients j\ and l\ one can, in principle, solve these 
two equations for ji and ^i. In practice, however, lack 
of precision in measuring the amplitudes Ap, Ai-a and 
Si-a is likely to produce large uncertainties. 



D. Yang- Yang anomaly 

The Yang- Yang relation (1.1) in the two-phase region 
(T < Tc) can be usefully rewritten as [2] 



Cv{T, p) = C*?* = {v/v,)Cp{T) + C^{T), (3.36) 
where v = V/N = 1/p and 

Cp{T) = v,T{d^p,/dT% C^{T) = -T{d^ii^/dT^). 

(3.37) 

The results (3.11) and (3.15) then yield 
Cp{T) = Ap\t\-" + B„ 



-'p 

-a-/3 



~ap\t\'^-" + ■ ■ 



+ hp\t\ 



h-a-l3+e 



C^{T) = A^\t\-'^ + B^ + ~a^\t\'-" + 



+ ci\tt 



-a-13-0 



(3.38) 

(3.39) 



where the various coefficients follow straightforwardly 
from (3.12)-(3.14), (3.16) and (3.17): see also K(3.103). 
Following Ref. [2] it is reasonable to define the strength 
of the Yang- Yang anomaly via 



lim — 



Cp + 



Ap -\- A^ 



By (3.32) this leads immediately to 

7^M = -i2/(i-i2). 



(3.40) 



(3.41) 



Note that 7?.^ depends only on _j2 (not on ji). Fisher and 
Orkoulas [2] estimated TZ^ for propane from experimen- 
tal data on the two-phase heat capacity and obtained 
TZfj, ~ 0.56: this suggests j2 —1.27. They also an- 
alyzed the heat capacity data for CO2 and estimated 
TZ^ ~ —0.4 (±0.3) which implies that j2 should be posi- 
tive but small (less than 1). On the other hand, simula- 
tions of the hard-core square-well fluid [6] indicate that 
TZfj, is small but negative, close to zero: correspondingly, 
j2 should be small but positive. 



IV. SPECIAL CRITICAL LOCI 

In this section we use the scaling formulation to obtain 
asymptotic expressions for various interesting critical loci 
that lie in the one-phase regions of the phase diagram. 
For convenience, we introduce a superscript index b de- 
fined so that: t = i for the locus ji = jic, say, the critical 
(3.35) isokyme; i = ii for the critical isobar, p = p^; i = iii for 



the critical isotherm, T ■ 
isochore, p = pc- 



Tc, and t 



iv for the critical 



A. Critical isokyme, p, = fic 

On the critical isokyme p = Pc (or p 
fields (1.4)-(1.6) reduce to 



p = p-kot- mop 
h = - j2P - kit 
i=t- jip 



rot - n3pt + 



m2p'^ 



r2t - n^pt + ■ 



mip^ 



■ Tit — riipt + 



0) the scaling 

(4.1) 
(4.2) 
(4.3) 
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Uh/\t 



-{hp 



kit) /\t — jip\^ now diverges, in general, when the critical 



Notice that the scaling variable y 
-jiP\ 

point is approached {i.e., when p, t — > 0) since A > 1 
for fluids. Hence, we need the expansions of the scaling 
functions W^{y) for y — *• oo. On using (2.3), (2.4) and 
(2.12), we obtain 



p = kot + mop^ + rot^ + nspt + 



+Qw^{u\h\)^'-"y^ 

+ QW^'>U^{U\h\)^' 
+ ahQW^'>U5c{U\h\)^^ 



l + w1t{U\h\) 



-l/A 



+ • 



(4.4) 



where an = sgn{h). Now i and h can be expanded in 
terms of p and t using (4.2) and (4.3) and the resulting 
equation may then be solved iteratively for p as a function 
of t. After some algebra, this yields the critical isokyme 
in the {p, T) plane as 

p\t) = \p\T)-p,]/p,kBn 

= P[t + + [l ± 5i |i|/3/A 



±6^|t|(^-i)/^ + --- + 6!,|t|«/^±6i|f| 



(/3+0)/A 



•••+ 6- + 6- |f I + ...]+..., (4.5) 

where ± refers to f ^ while 

p\ = ko, p2='r'o+ k^mo + koUs, 

4 = QW^{U\ki+j2ko\)^^-"^^^, (4.6) 

and the amplitudes b\, ■■■ of the correction terms are 
given in K(3.111). 

Note that the leading singidar exponent is (2 — a)/ A = 
l + ^~^, where (5 = A//3 is the standard critical exponent 
characterizing the critical isotherm. In the case of the 
{d = 3)-dimensional Ising universality class, S ~ 4.8 so 
that 1 + ~ 1.21. This implies that the curvature 
of p^{t), the pressure on the critical isokyme, diverges 
as T ^ Tc. Also by convexity — see Sec. II. A — the 
leading singular amplitude 4 is positive regardless of 
the sign of t. On the other hand, it transpires that the 
amplitude h\ of the leading odd correction contains the 
signum factor dh so that its sign depends on h. Later we 
will see that the sign of h can be determined from the 
mixing coefficient k\. 

To obtain the critical isokyme in the (p, T) plane, we 
first substitute (4.5) into (4.2) and (4.3) and express h 
and t as functions of t. Using these results in the defi- 
nitions (2.15) then yields p and s in terms of t. Finally, 
from (2.18), we obtain the density on the critical isokyme, 

= /ic, as 

p\T) = pe{l ± B^\tfl^ + Bii\t\^f^/^ + 4|t|(i-«)/^ + 

• • • ± + • • • + Bl\4P+'>^y^ + •••}, 

(4.7) 

where ± refers to /i ^ 0, while the leading coefiicients are 



B' = {1+ 5-^){l - i2)4/|fci + i2fco|, 

4 = -32{2 -a + /3)(B')V(2 - a)(l - 32), (4.8) 

and the further coefficients are presented in K(3.113). 

Note that the leading singular exponent is/3/A = (5~-^ 
which is less than /3 so implying that the the critical 
isokyme in the (p, T) plane is significantly flatter than the 
coexistence curve. Clearly, therefore, the critical isokyme 
helow Tc lies outside the coexistence boundary, i.e., en- 
tirely in the one-phase region as is to be expected. 



B. Critical isobar 

At p = pc (or p = 0), the scaling fields (1.4)-(1.6) with 
(3.22) reduce to 



kot- fx - qofi^ - rot^ - vofit + ■ 



p, — kit — q^fi^ — r2t^ - + • • 



t = t — lifi — qijj? — rit^ 



Viflt - 



(4.9) 
(4.10) 
(4.11) 



As the critical point is approached along a general lo- 
cus {i.e., t, jj — > 0), the scaling variable y oc h/\i\^ 
(/U — kit)/\t — lifi\^ again diverges. Using the large-y 
expansion (2.12) for p and rearranging (4.9) yields 



jl= - kot- qojl^ - rot"' - voixt 



+ 



- QW^(?7|/i|)(2-")M 1 + w°ii{U\h\)-^/^ 

-0l^i^)C/4c([/|/l|)('-"+')/^[l + ---] 

- QW^^U5c{U\h\)^^-"+^'^/'^[l + ■■■]. (4.12) 

Solving this equation iteratively for fi in terms of t with 
the aid of (4.10) and (4.11), expresses the critical isobar 
in the (/x, T) plane as 

/i"(<) = [M"(T)-Mc]/fcBTe 

= -kot + jlit' + • • • + 4|t|(2-«)/^ [1 ± 

± 6^-|f|(^-i)/^ + • • • + }}^\tf^ ± }^l\t\''fi+'^'^ 
+ ■■■ + b^^\tf^'^ + fe^'ltlC^+^^'Z^ + •••], (4.13) 

where ± again refers to t ^ 0, while the leading coeffi- 
cients are 



A" = -ro + koVo - kgqo, 
A'; = -QW^{U\ki+ko\)^^-"^^^, 
6« = -(2-a)4VA(fci + fco), 



(4.14) 



and the further correction amplitudes are given in 
K(3.179). 

Note that the curvature of /i"(t) diverges at the critical 
point with the same exponent as does the critical isokyme 
in the {p,T) plane: see (4.5). The critical isobar also has 
the same sign of curvature above and below Tc in the 
{jj., T) plane. By thermodynamic convexity, the leading 
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singular amplitude A" is negative. The amplitude 5" of 
the leading odd correction term again changes its sign 
depending on h. 

In the (p, T) plane the critical isobar can be obtained 
via the same route used above for the critical isokyme. 
The result is 

p"{T) = pe{l ± + -Bj'ltl^^/^ + + 

■ • • ± + • • • + + ...}, 

(4.15) 

where ± here refers to /U ^ /Uc while the leading coeffi- 
cients are 

B" = -(1 + 5-i)(l - j2)A^/\h + ko\, 



[(2 - a)j2 + P]{B"r/{2 - a){l - h), (4.16) 



where the further coefficients are to be fomid in K(3.121). 

Notice that the leading singular behavior matches that 
of the density on the critical isokyme as given in (4.7). 
For small j2 (< 1), the leading amplitude must, by 
convexity, again be positive. The ratio between in 
(4.7) and B" is simply 



BVS" = |(A;i+j2A;o)/(fci + fco)l 



C. Critical isotherm 



0/A 



(4.17) 



When T = Tc so that t 
reduce to yield, now. 



0, the scaling fields again 



p = p- jl- mop^ - goA^ - "oPM + 



(4.18) 



and similarly for h and t. Once more, the scaling variable 
y diverges, in general, on approach to criticality. Using 
(4.18) and the large-y expansion (2.12) for the scaling 
functions yields an equation for p in powers of jl, h, t 
and p which, as before, can be solved iteratively with the 
aid of the reduced expansions for h and i to obtain jo as a 
function of fi. The result for the pressure on the critical 
isotherm is 

r'(/i) = [p'"(M)-Pc]/PcfcBTc 

+ Af\jl\^'^-"'y^ [l ± bf\fif'^ ± • • •] , (4.19) 
where ± refers to /i /ic, while the coefficients are 



p2 =mo + go + «o, 

6i" = -j2(2-a)4"/A(l-j2), 



(4.20) 



and the spectrum of higher order terms matches that in 
(4.13): see also K(3. 127). The leading singular exponent 



is again (2 — a)/A = 1 + 5~'^ , implying that the critical 
isotherm in the (p, ji) plane has a divergent curvature at 
the critical point, and thermodynamic convexity ensures 
4" > 0. 

To obtain the critical isotherm in the {p,p) plane, we 
first invert (4.19) to obtain jl k p - and 
then use this to express the generalized densities, p and 
s, as functions of p. Using (2.18) finally yields the critical 
isotherm in the form 

f^\p) = Pc {l ± E^'W"^ + Bflpl^A + ^iii|^|(l-a)/A 

+...±Bf\p\^^+'y'^ + -.- 

+ + •••}, (4.21) 

where ± again refers to /i ^ /ic while 

= (1 + 5-1)1(1 - ji)/(l - j2)|''/^4", (4.22) 

and Bf oc -{B"'f: see K(3.131). As expected for the 
critical isotherm the leading exponent is again /3/A = 
5~^. Assuming that the pressure- mixing coefficient j2 is 
small (< 1), we find B'" > owing to convexity. Now 
when j2 — Q, the scaling field h can be approximated by 
p which implies that the critical isotherm in the {p,p) 
plane approaches the critical point from higher density 
above p^-, while from lower density below p^- this accords 
with the observed standard behavior. 

D. Critical isochore 

Below Tc in normal fluids the critical isochore in the 
(/U, T) and {p, T) planes coincides with the phase bound- 
ary, HrriT) and PaiT) (but see [23] for exceptions in cer- 
tain models); however, the behavior above is of general 
interest. When p = pc, the result (2.18) for the density 
leads to 

= (1 - j2)QU\if [Wl'{y) + U4c\i\'wi^^'{y) + ■■■ 

+ j2{j2-l)Q^U^\i\^^[wl'{y) + ---f 

+ {h +ji)Q\i\'-'^[AyWl'{y) - (2 - a)Wl{y)\ 
+ ■■■, 

(4.23) 

where the primes denote differentiation with respect to 
y. We need to solve this equation for y as a function of t; 
but on the critical isochore we expect y ^ when t ^ 0. 
Hence we should now use the small-y expansions (2.8) 

for the scaling functions W!l{y), w!i^\y), etc. The re- 
sulting equation in powers of y may be solved iteratively 
to obtain 



with coefficients 



(4.24) 
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Fi = 



{2-a){h+,h)Wl 



+0 



2{l-j2)Wl^U 



The scaling fields p and h along the critical isochore 
can now be found as follows: first, on combining (2.3), 
(2.4) and (2.8) with (4.24), we find 



;5 = 0|f|(2-") [W% + Wl,Y,'\t\ 



2-2Q-2/3 
.(5)i 



+ 



+ •••]; 



(4.26) 



then, from the definition (2.1) of the scaling variable y 
we get 



/l=(Fl/C/)|tT-«+T [l + y2\if + . 



(4.27) 



finally, by expressing the scaling fields in terms of p, fi 
and t via (1.4)-(1.6), wc can solve these two equations 
iteratively for p and fi as functions of t. After some al- 
gebra, this yields the critical isochore in the {p, T) plane 
as 

f-{t) = \f-{T)-p,]/p,k^T, 

= Pa,lt + p,,2t' + • • • + [l + C^^X 

+ 4;iir"" + • • • + + • • • ] 

+ B7|i|^-"+T[l + b^;\t\^ + ■■■], (4.28) 
where the leading amplitudes are 



.\1-a 



+01 



B'; = Yx\t\^-''^^ lU, (4.29) 



while Pcr,i) Per, 2 5 T and the correction amplitudes are give 
in (3.12), (3.14) and K(3.139). The amplitude J^^^^ is 
positive by convexity. 

In the T) plane, the critical isochore is given by the 
closely analogous form 



A'"(f) = [m'"(t) - McI/^bTc 

= /i<T,ii + Aa,2t^ + ---+i2^jr|i 



ivu|2-a|^^^iv|^|e 



■a2;iir-" + --' 



4;i*r""-^+'=+---] 



(i+j2)i37ii|i-"+^[i+6^-itr 



(4.30) 



in which /x^.i and p,c!,2 are given in (3.16) and K(3.78), 
while the coefficients A^^ , a5^, etc. are the same as in 
(4.28). 

When the pressure mixing coefficient ji vanishes, the 
leading singular term here vanishes; but in that 

case the third-derivative of /Lt''^(T), the chemical potential 
on the critical isochore above Tc, diverges at criticality, 
since one typically has 2<1 — q; + 7<3. In view 
of these results one might also define a magnitude for 
a Yang- Yang-type of anomaly by comparing the leading 
singularities in p and \i on the critical isochore above Tc', 
but this naturally leads to the identical ratio! One may 
thus write 11+ = -ji/il-ji) = T^ix- see (3.41). 



E. fe-susceptibility loci 

In this section we focus on the novel fc-loci defined 

in the one-phase region via the isothermal maxima of 



X 



(fe) 



Xl p in the (p, T) plane [6], where x = {dp/d^) 



T- 



see also (2.22). If one considers x^*'^ as a function of /U 
and T, the maxima of x^'^^ at fixed T satisfy 



V J J, 



dp 



0. (4.31) 



This leads to the condition 



or 



see (2.14) and (2.21). By using (2.29), one can obtain 
the /t-derivative of Xivjv and express the fc-locus in terms 
of the generalized susceptibilities introduced in (2.28). 

After some algebra, one finds 



Xhhh + (1 + j2)pXhhh - eo{k)Xhh - ^(^iXhht 



where, by convention, Xhhh 

{d^p/ dli^dt), while 



•• = 0, 
(4.33) 

{d^p/dh% and Xhht = 



eo(fc) = 3j2 + fc(l - 32) and 64 = {h 



■Ji)/(1-J2). 

(4.34) 



The condition may now be converted to scaling form by 
using (2.3) and (2.4); then by employing the small-y ex- 
pansions (2.8) one can solve to obtain 



y = Uhl\tf = eo{k)Y^\tf + Y2\t\ 



\-a-l3 



+ ■ 



(4.35) 



which might be compared with (4.24) for the critical iso- 
chore. The coefficients arc given by 



Y, = \QU{Wl,f/Wl„ 



?2 = lie^WlJUWl^, 

(4.36) 



and wc may note, for its future significance, that the 
leading term in (4.35) vanishes identically for the special 
k value 



fcopt = -3j2/(l -32) 



(4.37) 



At this point the scaling field p can be expanded in 
powers of y and \t\ via (2.3) and (2.8) and then, via (4.35), 
wholly in terms of Finally, by using (1.4)-(1.5), and 
rewriting (4.35) as an expansion for h one can solve for 
p iteratively as a function of t. After some algebra, we 
find 



(i) = Pa,it + p^^^t" + • • • + 4^=) 



P 



+ ■ 



(4.38) 
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where the leading amplitudes are 

= QW%\rf-^ [1 + eo{k)iW%f/6W%W'_ 



1 

+4j : 



(4.39) 
(4.40) 



while Pcr,i, Pa.2, T and the remaining coefficients are given 
in (3.12), (3.i4) and K(3.153). 

The fc-loci in the (/i, T) plane can now be obtained by 
substituting this result into the scaling fields and using 
(4.35) once more. The result is 



/iW(f)=/i<,,lt + /i.,2t' + ---+A«|t|' 



+ 4!)|ir + ... + a.r/|i 



(fc) 1^1/3+65 



3/i I 



(4.41) 



where /ig-.i and /io-,2 are given in (3.16) and K(3.78) while 
the principal amplitudes are 



A^^^ =j2A<^^+eoim\T\^-yU, 
i?« = (l+j2)B('=\ 



(4.42) 



and ag) - .724'^«SV4'^ for z = 1, 2, 3. 

For practical purposes the form of the A;-loci in the 
density- temperature plane is of most interest. To that 

end, note that the generalized densities, p and s, can be 
written using the above results as 

p = 2Qu[eo{k)YiWl^\i\^'' + Y^Wl^lil^-" + ■■■ 

+ eo(fc)ri[/4cH^|t^|t?''+' + --- 

+ ^U5cW'f^\tf+'' +■■■], (4.43) 
§ = Q|iT-« [(2 - a)W% - ^[eoik)Y,]'Wl,\i\''^ + •••], 

(4.44) 

where t k, rt with t defined in (3.14). The fc-loci in the 
[p, T) plane are given by 

pi^) (i) = 1 + (1 - j^yp + j2(j2 _ i)p2 _ + + . . . ^ 

(4.45) 

so that finally the fc-locus varies as 

p«(r) = pc {i + Bf\t\^^ [1 + h,\tf + •••]+ Af\t\'-'^ 

+ 4'=)* + • • • + + •••}, (4.46) 

where the coefficients are 



2/3 



Bf^ = 2eo{k){l - 32)Y^QUWI^\t 

Af ^-{h+n)Q\T\'-'^ 

x[{2-a)W% + Hwl,f/W%] 
= (l-i2)QC/C/5cW^|'?|rr^^ 



(4.47) 

(4.48) 
(4.49) 



while the expression for A\ ' is rather complicated. What 

is significant is that the leading amplitude, -B^'^^ varies 
linearly with k and vanishes identically when k takes the 
"optimal value" fcopi = 37^^ given in (4.37), while the 
other coefficients exhibited do not vary with k [despite 
the superscript label (A;)]. For k = fcopt we may say, 
loosely, that the fc-locus points most directly to the crit- 
ical density pc- Indeed, for this reason examining the 
fc-loci may be of value in analyzing both experimental 
and simulation data. Furthermore the relation to the 
Yang- Yang anomaly ratio is again revealing and sugges- 
tive. 

Orkoulas, Fisher and Panagiotopoulos [6] examined 
the fc-loci for the hard-core square-well fluid using grand 
canonical Monte Carlo simulations. They observed that 
the fc-loci for different system sizes settle down and be- 
come independent of size at high enough temperatures. 
Within the precision attainable these loci can be consid- 
ered as the true fc-loci (for the thermodynamic limit). 
However, when T T^, the finite-size loci clearly devi- 
ate from the limiting behavior. For the data in hand the 
finite-size effects become evident when t = {T — Tc)/Tc < 
0.1. 

To estimate the fc-loci near the critical point, we have 
fitted the data for t > 0.1 with the formula (4.46) retain- 
ing only the coefficients b'^'^ , A^'^-' and A^"^ , while adopt- 
ing Ising values for the exponents and taking pc = 0.3067 
and Tc = 1.2179 [6]. Some of these estimates are pre- 
sented in Fig. 1. Similarly, some of the estimated fc- 
susceptibility loci for the restricted primitive model [25] 
are shown in Fig. 2. For comparison, Fig. 3 presents the 
fc-loci for the van der Waals fluid: see Sec. V. Note that 



0.25 



t 




-0.25 



FIG. 3. Selected fc-susceptibility loci for the van der Waals 
equation. The thick solid line represents the coexistence curve 
while the thick dashed line is the diameter. The (fc = 0) locus 

is the vertical solid line, while the fc = 1, |, i and fc = — 1 



loci are portrayed by long dsished, short dashed, dotted, and 
dot-dashed lines, respectively. Note that the fc = f locus has 
the same slope at criticality as the coexistence curve diameter. 

for the van der Waals fluid all the fc-loci approach the 
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critical point linearly, which, of course, is consistent with 
the classical exponent equalities 2/3 = 1 — a = 1. 



F. fc-heat-capacity loci 



It is also interesting to examine the fc-hcat-capacity loci 
or Cy ■'-loci defined by points of maxima of the modified 



specific heat 



Cf\T,p)^Cv{T,p)/p^ 



(4.50) 



in the (p, T) plane, where Cy is the constant-volume spe- 
cific heat. These are clearly quite analogous to the k- 
susccptibility loci discussed in the previous section; but 
they have not, as yet, been used in any simulations. The 
C^^-loci can, in principle, be obtained in a way similar 
to that used for the fc-susccptibility loci by starting with 
the relation (2.27). However, when one takes a deriva- 
tive of (2.27) with respect to n at fixed T, the expression 
becomes complicated and difficult to handle. Therefore 
we outline a different, canonical approach. 

The required maximal points in the {p,T) plane at 
fixed T satisfy 



kCv +p{dCv/dp)j, = 0. 



(4.51) 



To find a convenient expression for Cy (T, p) , we consider 
the Helmholtz free energy density f{p,T) = pp{p. T) — 
p{p, T). In terms of the reduced variables p, p, and p one 
has 

f{p,t) = [f-u] /p.kBn 

= {pc/kBTc)Ap + pp - p with Ap = p-1. 

(4.52) 

The reduced specific heat, (2.25), is then 

Cvip, t) = {pn/kBPcT)Cv = - (d^f/dt^)^ , (4.53) 

while the fc-locus equation, (4.51), becomes 

-{k + l)Cv+p{dCv/dp)^ = 0. (4.54) 

To solve this for p as a function of t, we first expand 
/, about the critical density in powers of Ap for < > 0. 
If the expansion coeflacients are fo{t), fi{t), ■ ■ ■, we can 
rewrite the locus equation as 

-{k + 1)^' + fi' + [{k - l)fi' + 2f^]Ap + 0{{Apf) = 0. 

(4.55) 

To expand the coefficients fo(t), etc., in powers of t, 
notice first that, from (4.52) with p = 1 or Ap = 0, we 

have 



(4.56) 



where /i'^(t) and p^'^{t) represent the variation of p and p 
on the critical isochore as obtained in (4.28) and (4.30). 
Then the relation p{T,p) = {df/dp)T yields 

Mt) = {fic/kBn) + r{t). (4.57) 

Finally, we have 

Mt) = i {dp/dp\\^^^^ = ^XNNit-^P^ Pc). (4.58) 

To obtain the reduced susceptibility, Xnn, on the critical 
isochore, we may use (2.29) and the previous results in 
Sec. IV.D. After some algebra we obtain 



XiVJV = 2(l-i2)WW?2N" 



i + Ui,iwi^^/wl^)\Ttf + - 



(4.59) 



where r was defined in (3.14). Note that the pressure- 
mixing coefficient j2 first enters in a t^~" correction; but 
that is of higher order than the term retained here. 
On taking the reciprocal and differentiating twice with 
respect to t, we finally have 



2f^{t) = Df-^ [1 - det^ + ■ 
where the coefficients are 



(4.60) 



D = 



7(7-l)|Tp 



2(l-i2)WH^^2' 



7(7-l)W^^2 



(4.61) 



We are now in a position to solve (4.55) iteratively for 
Ap as a function of t by using (4.28), (4.30) and (4.60). 

One finally obtains the /c-heat-capacity or Cy -locus in 
the form 



p'-c\t) = P. 



1 



(4.62) 



where, recalling (4.29), (3.12) and K(3.68)-(3.78), the 
amplitudes are 



(fe) 

c 

(k) 



(2-a)(l-a)[l + fc(l-j2)]47A 

[Pa,2 + k{p„,2 - Pa.2)]/D. 



(4.63) 



Note first that the exponents 2/3 and 2/3 -fa are numer- 
ically very close (since a ~ 0.1), so the two terms derived 
compete strongly near the critical point; furthermore, a 
Ifl"^"" term also appears in the full expression. However, 

the leading amplitude ^ does vanish when 



k = kc 



(4.64) 



This value can be compared to 37?.^ for the A;-loci: see 
(4.37). However, one must keep in mind that it will be 

more difficult to resolve the 'optimal' k value here, com- 
pared to the fc-susceptibility loci, since the term varying 
as t'^l^+°' does not vanish at k = kc- 
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FIG. 4. The fc-heat-capacity loci for the hard-core 
square-well fluid in a finite periodic cube of side L derived 
from previous simulation data [6]. The dotted curves are the 
loci for L* = L/cr = 6, where a is the diameter of the hard 
spheres; the dashed lines are for L* = 9, while the solid lines 
arc for L* = 12. Note that since the systems are finite these 
loci (and the fc-susceptiblity loci) extend below Tq. 



In order to gain some impression of the behavior of 
these Cy ^-loci, we present in Fig. 4, some of the fc-heat- 
capacity loci obtained for the hard-core square-well fluid 
in finite systems. Note that since the exponents 2/3, 2/3-1- 
a, 1 — a, etc., arc closely spaced, it is not feasible to 
extract reliable estimates of the thermodynamic limiting 
loci from such finite-size data. 



V. CRITICAL LOCI IN CLASSICAL THEORY 

In this section, as concrete, albeit rather special, ex- 
amples of the various critical loci discussed above, we 
consider the classical theory of a liquid-gas critical point. 
In particular, the van der Waals equation will be used to 
provide quantitative illustrations: it reads 



p = pkBT/{l - bp) - Ap'' 



(5.1) 



where h and A are constants which measure the molec- 
ular size and the strength of the attractive interactions, 
respectively. As well known, one has 



/cbI; = 8A/276, pc = 1/36, Pc = A/271? 



(5.2) 



Less well known is the behavior when T ^ of the liquid 
and vapor densities and of the phase boundary /x^ {T) — > 
-A/b: see K(App. D). 



A. Phase boundaries 

The Helmholtz free energy (or its appropriate analog) 
in a classical or Landau theory may be taken as analytic 



throughout the critical region. Accordingly, to formulate 
the theory we expand the free energy density around the 
critical point in terms of the order parameter 

m={p- pc)/Pc (5.3) 

and the temperature deviation t — {T — Tc)/Tc, as 



/(T,p) = 5^5^a,•fcm^^^ 

j=0 k=0 

with = ^30 = 0, > 0, (5.4) 

where the conditions stated serve merely to ensure nor- 
mal critical behavior. The explicit values of the leading 
coefficients ajk for the van der Waals equation (5.1) are 
derived in K(App. C) where it is seen that a special fea- 
ture is that all the cubic coefficients, as^k {k = 0, 1, 2, • • •), 
vanish identically, as do many higher order coefficients 
such as aj2, aj^, etc., for all j > 2. Of course, these 
features should not be expected to hold in real systems 
or in more realistic models even when a classical descrip- 
tion of criticality may be warranted. The leading nonva- 
nishing van der Waals coefiicients are reproduced in the 
Appendix here. 

Prom (5.4), the chemical potential and the pressure 
can be expanded using 



p{T, p) = {df/dp)^ , p{T, p)=ptx- f. 



(5.5) 



In the two-phase region, the liquid and vapor phases, 

with densities puq = pdl + mnq) and pvap = Pc(l+?72vap), 
respectively, must have the same chemical potential, 
Pa{T), and pressure, Pa{T). Solving these two condi- 
tions for piiq and pvap by using (5.3) and (5.4) yields the 
desired phase boundaries. The detailed calculations are 
presented in K [26]; the results are 



m(T) 
mo(T) 

lia{T) 



^(Wliq • 
|(miiq 



m 



vap; 



Alt + Ait^ + 0{t^), (5.6) 
[l + Cit + 0{t'')] , 
(5.7) 

: + Pit + p2t'^ + Pst"^ + 0{t'^), (5.8) 
■■Pc+Plt + P2t^ +P3t^ + 0{t^), (5.9) 



where the amplitudes Ai, A2, B, etc. are expressed in 
terms of the coefficients ajk in the Appendix. The re- 
sults agree with those of Sengers and Sengers [13] who, 
however, give them only to one order lower in t. 



B. Critical isochore and analytically continued loci 

In addition to examining (i) the critical isochore above 
Tc one may, for critical points describable by classical the- 
ory, always consider the analytical continuation toT > 
of various loci otherwise defined only for T < T^. Specif- 
ically, we will study: (ii) the continued coexistence curve 
diameter, p{T); (iii) the continued saturation chemical 
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potential, ^a{T)] and (iv) the continued vapor pressure 
line, pu{T). For nonclassical critical behavior, the analy- 
sis of the previous section demonstrates that, in general, 
these last three continuations cannot be uniquely defined 
since all carry singularities at Tc. 

We address first the appearance of these loci in the 
(p, T) plane: see Fig. 5. The critical isochore, (i), is triv- 
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FIG. 5. Various critical loci for the van der Waals cquar- 
tion of state in the (p, T) plane, where m = (p — pc)/pc- 
The coexistence curve is drawn with a thick solid line, (i) 
Critical isochore; (ii) coexistence curve diameter and its an- 
alytic continuation into the one-phase region; (iii) analytic 
continuation of fj,a(T); (iv) analytic continuation of Pa{T); 
(v) (fe = l)-susceptibility locus. 

ial; the continued diameter, say Pii{T) [ = p{T)], follows 
directly from (5.7) and has a critical slope 



{dpii/dt)c = Aipc, 



(5.10) 



which is negative for the van der Waals equation: see Fig. 
5 for which the quadratic terms in t were also computed. 

To determine Piii(T), the density locus on which the 
chemical potential is the analytic continuation of the 
phase boundary, p„{T), we may substitute (5.8) in the 
full expansion for p{T,p) that follows from (5.4) and 
(5.5). This gives an equation connecting m and t which 
is easily solved for m in powers of t although the p2t^ 
term in (5.8) is needed even in linear order. One finds 

{dpiii/dt)c = Pc (a2ia5o - 2a3ia4o) /8a|o- (5-11) 

The locus piii(r) can be regarded as an effective line of 
symmetry along which the chemical potential is analytic. 
However, as seen in Fig. 5, this locus differs from the 
analytic continuation of the coexistence curve diameter, 
Pii{T) — another natural candidate — even in the lowest 
order in t. 

The locus piv{T), along which the pressure is the con- 
tinuation of the vapor pressure curve, p^{T), can be 
found in an analogous way. Substituting (5.9) into the 
expansion of p following from (5.5) leads to 



(dpiv/rft)c = Pc (01210140 + 0121050 - 2a3ia4o) /8a|o, (5.12) 

which, in fact, differs from both (5.10) and (5.11): see 
Fig. 5. 

In the {p,T) plane one can find results for the same 
four loci by using (5.4) and the expressions found for the 
loci in the (p, T) plane. In all cases the initial slopes are 
the same, that is {dpJdT)^ = pi for i = i - iv: sec Fig. 
6. However, the initial curvatures differ, being given by 




FIG. 6. Critical loci for the van der Waals equation in the 
(p,T) plane: pcAp/pc is plotted with, for illustrative pur- 
poses, Ap = p — Pc — (pi -|- |)t [see (A5) and (A13) for pc 
and pi]. The graphs include only the second and third order 
terms in t. The phase boundary, pa(T), is drawn with a thick 
line. The labeling (i)-(v) is the same as in Fig. 5. 

PcPi'c = 2ai2 and 

PcMi'i.c = 2ai2 - a2i(2a3ia4o - 2a2ia5o)/2a4Q, (5.13) 
PcpS'ii.c = 2ai2 - 021(2031040 - 021050)72040, (5.14) 

2oi2 — 021(2031040 — O21O50 — O2l04o)/204Q, 

(5.15) 



PcPiv.c 



where the primes denote differentiation with respect to 
t. 

Similarly, the loci in the (p, T) plane can be obtained: 
see Fig. 7. Again, all have the same initial slopes while 
the curvatures are distinct as follows from 



p-'c = 2ai2 - 2oo2 
p"i c = 2ai2 - 2oo2 



40' 



Pm.c = 2ai2 - 2oo2 



Piv,c = 2ai2 



021(2031040 - 2a2ioi5o)/2o 

O2l(2o3iO40 — a2l05o)/204o, 
2oo2 — 021(2031040 — 021040 
-O2l05o)/204Q. 



(5.16) 
(5.17) 
(5.18) 

(5.19) 



Notice that in Figs. 6 and 7 the critical values as well 
as conveniently chosen terms linear in t have been sub- 
tracted from each locus. Thus one can clearly resolve the 
differences in curvature and observe that the sequence of 
loci, from top to bottom, is the same as in Fig. 5. 
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FIG. 7. Critical loci for the van der Waals equation in the 
(p, T) plane, where Ap =p— — (pi — 2)f [see (A8), (A5) and 
(A12)] and the graphs are correct only up to third order in 
t. The vapor pressure curve, Pct{T), is represented by a thick 
line. The labels (i)-(v) have the same meaning as in Figs. 5 
and 6. 



C. fe-susceptibility loci 

Since the density increases monotonically with fi at 
fixed T, the maximal condition (4.31) specifying the k- 
susceptibiHty loci can be rewritten as 



(k) 



dp 



= 0, 



(5.20) 



and, thence, in terms of the free energy density f{T,p), 
in the simpler form 



(5.21) 



Substituting the expansion (5.4) and solving for m iter- 
atively yields the fc-loci in the (p, T) plane generally as 

pW(T) [I + m,{k)t + m2ik)t'' + 0{t'')] , (5.22) 

with coefficients polynomial in k, namely, 

mi(fc) = -(fca2i + 3a3i)/12a4o, (5.23) 
m2{k) = -yj {ka22 + 3a32 + 3[{k + 1)031 + 4a4i]mi(fc) 
+ 6[{k + 2)040 + 5a5o][mi(k)f} /ma. (5.24) 

For the van der Waals equation these results yield 

mi{k) = -lk, m2{k) = -lk{k^ + k - Z) (5.25) 

(see Appendix). Note that when fc = 0, the coefficients 
mi and TO2 both vanish; furthermore, an exact calcula- 
tion shows that all the expansion coefficients vanish iden- 
tically at fc = so that one has p'^^^{T) = p^. see Fig. 3. 
We also find that the slope of the fc-locus at criticality 
becomes equal to that of the coexistence curve diameter 



when fc = |. In Fig. 3 the fc-loci for the van der Waals 
fluid at several values of k have already been presented. 
Comparison with Figs. 1 and 2 reveals that the behavior 
of these loci is closer to those for the hard-core square- 
well fluid than for the restricted primitive model. Figs. 
5-7 show how the (fc=l)-susceptibility locus appears in 
the {p, T) and {p, T) planes and relates to the other van 
der Waals loci: see plots labeled v. 



D. Yang- Yang relation and some extensions 

It is natural to ask how, if at all, a Yang- Yang anomaly 
might appear in a classical theory. To that end we dis- 
cuss, in this section, the Yang- Yang relation (1.1) on gen- 
eral linear loci in the (p, T) plane and also derive an anal- 
ogous relation for isotherms. 



1. On the critical isochore 



On the critical isochore p = pcin the two-phase region 
below Tc, it is convenient here to define the functions 



M-{T) = -p,T- 



(5.26) 



where the factors have been introduced simply for 
dimensional convenience. [Compare with (3.36)-(3.37).] 
Above Tc the functions C+{T), V+{T) and M+{T) may 
be defined on the critical isochore in precisely the same 
way. The Yang- Yang relation then implies = + 
M^. In classical theory we may expand in powers of t 
for f ^ to obtain 



(5.27) 



and likewise for V^iT) and M.'^iT). The leading ampli- 
tudes above Tc are simply 



-2a, 



■02, 



7'+ = 2(ai2-ao2), Ml 



-2ai2, 
(5.28) 



while the amplitudes of the terms linear in t are 

Cr =6(P3 -PcMs), Ci+ = -6ao3, :Pr=6p3, (5.29) 
■p^ = 6(ai3 - aos), = -QpcPs., = -6013, 

where pz and pz are given in the Appendix from which 
one readily sees that the slopes are discontinuous across 

Indeed, it is well known that the specific heat Cy(T) 
on the critical isochore exhibits a finite jump at a classical 
critical point which, in fact, is of magnitude given by 



PcTcACy = AC = C- 



■C+- 



a^i/a4o. (5.30) 
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It transpires, however, that in general both the second 
derivative of the pressure and of the chemical potential 
have a finite discontinuity on the critical isochore. This 
arises from the presence of the odd coefficients 031 and 
050 (even though 031 vanishes "by accident" for the van 
der Waals equation) as follows from 

AM = M~ - A1+ = 5021(2031040 - a2ia5o)/alo, (5.31) 

with, of course, AV = AC - AM. 

These results are illustrated for a van der Waals fiuid in 
Fig. 8. Clearly, the variation of both the pressure and the 
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FIG. 8. The specific heat Cy(T, pc) (solid curve) on the 
critical isochore of a van der Waals fluid (times i in units of 
/sb), compared with the corresponding contributions, Cp{T) 
and Cfj.{T) (dot-dashed and dashed plots), due to the iso- 
choric variation of pressure and chemical potential: see (3.37) 
and (5.26). Note that the standard kinetic contribution to 
the total specific heat, namely ^dk^ (in d dimensions), arises 
entirely from — — T/i"; the plots beneath Tc are correct 
only to leading order in t. 

chemical potential contribute to the discontinuity in the 
specific heat on the critical isochore. The relative degree 
of the two contributions may be gauged by evaluating 
the ratio 

AM ... ^^_3^^ 



AC 



^021 



050 

040 ' 



which, as the notation suggests, might be regarded as an 
effective Yang- Yang ratio for classical systems. Indeed, 
if TZfj, vanishes, as in the case of simple lattice gas mean- 
field models, the isochoric second temperature derivative 
of the chemical potential becomes continuous through the 
critical point: then, the only contribution to the discon- 
tinuity in the specific heat on the critical isochore arises 
from the variation of the pressure. 



2. On linear density loci 

In light of the various critical loci that approach the 
critical point linearly in the {p, T) plane of classical sys- 



tems (see Fig. 5), it is rather natural to extend the Yang- 
Yang relation to the general linear locus 



p = Pe(l+rt). 



(5.33) 



To that end consider, first, the second temperature 
derivatives of the pressure and the chemical potential 
near criticality along this locus. For T < Tc, there is 
no dependence on r because, since the exponent /? is less 
than unity, the linear locus (5.33) always lies in the two- 
phase region when T Tc — . Hence the results (5.29)- 
(5.31) still apply. In the single-phase region above T^, the 
previous definitions may naturally be extended by taking 



Mf{T) 



(5.34) 



From (5.4) and (5.5) these functions may be expanded 
straightforwardly to yield the limiting values 



(5.35) 



P+c = 2(2ai2r + 012 - O02), 
7W+c = -2(012 + 2o2ir). 

In an obvious notation, we then obtain 

AVr = 021(021040 - 2031O40 -I- O2i05o)/2o|o " 4a2ir, 

(5.36) 



AMr = 4:a2ir + 021(2031049 - 021050)7204 



(5.37) 



When r = 0, i.e. along the critical isochore, these results 
of course agree with (5.31). 

It is now interesting to define r-p as specifying that 
locus along which the pressure has a continuous second 
derivative at criticality and similarly for r^. From (5.36) 
and (5.37), we thus find 

rM = (021O50 - 2o3i04o)/8o4o = r-p - §021/040. (5.38) 

For the van der Waals equation these expressions yield 
r-p = 0.8 and rj^ = —0.2. On the other hand, from 
(5.10), the locus of the analytically continued coexistence 
diameter is specified by 



(021O50 - 031040)74040, 



(5.39) 



which takes the value f = —0.4 for the van der Waals 
equation. Evidently, all three loci are distinct! Note 
especially, however, that rx is equal to the slope of the 
effective line of symmetry at the critical point: sec (5.11). 
Hence, as already shown by MulhoUand [28], this puta- 
tive line of symmetry for the van der Waals equation 
differs from the analytic continuation of the coexistence 
curve diameter; but it is, rather, the locus on which the 
efi^ective Yang- Yang ratio, 7t^, vanishes! 

Now, extending the Yang- Yang relation to the general 
locus (5.33) leads to 



pc^^ py 



T^P^Kt 



T-^[Vt{T)+Mf{T)\, (5.40) 
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where p = Pc(l + ^t) is understood and, as before, Kt 
denotes the isothermal comprcssibihty: see K(App. F). 

In classical theory, this relationship may be verified 
straightforwardly using (5.4), and subsequently derived 
expressions for pCy/T and 1/ p^Kx, on (5.33). Explic- 
itly, one finds the initial derivatives of (T) and (T) 
at criticality to be cubic polynomials, namely, 

"P+i = 6 [ai3 - ao3 + 2a22r + (021 + 3a3i)r^ + 4a4or^] , 

(5.41) 

M+-^ = - [6ai3i + 2(ai2 + 6a22)r + (4a2i + ISas^r^ 



(5.42) 



which, of course, satisfy the generalized relation (5.40) 
with 



pl/p^KT^2a2it + 0{t'). 



(5.43) 



3. Yang-Yang-type relation on the critical isotherm 

An important locus in the (p, T) plane is the critical 

isotherm, T = Tc- In analogy to the Yang- Yang relation 
for isochoric variations one can derive a corresponding 
isothermal relation, namely. 



(5.44) 



see K(2.97). 

To apply this to the critical isotherm, recall that the 
generally expected critical behavior is 



l/plKT^D^\Apf-\ for Ap = p-p,^0±. 

(5.45) 

In classical theory one has 5 — ?> and finds D+ = D~ = 
12a4o/Pc- More generally one might define the pressure 
and chemical potential contributions in (5.44) via 



v,{p) = (aVV)r , M,{p) = - {d^p/dp'). 

at T = Tc. 



(5.46) 



It is then natural to ask, as on the critical isochore, how 
the two contributions — one from p, and one from p — 
contribute to the overall singularity in (5.45) and whether 
that might throw any further light on the pressure- mixing 
coefficients ji or 72- It turns out, however, that matters 
are very different! In fact, by appropriate integration of 
(5.45), one can obtain the leading singular behavior of 
both p and p on the critical isotherm. Then one finds 
that the leading singularities of Vpip) and A4p{p) both 
vary as |Ap|*~^ and are, thus more singular tha,n is their 
sum! Indeed, these leading terms must cancel exactly in 
(5.44) so that the behavior (5.45) for 1/Kt appears only 
as a net correction term. Needless to say, the classical 



theory fits in with this description although, since S is 
an odd integer, there is ambiguity in defining appropri- 
ate "singular contributions" . More generally, while there 
may well be effective ways in which an estimate of the 
degree of pressure mixing can be found from isothermal 
observations, we have not identified a good candidate. 



VI. SUMMARY 

We have carefully formulated a full or "complete scal- 
ing theory" for asymmetric fiuid criticality that, in par- 
ticular, incorporates pressure mixing in the basic linear 
and nonlinear scaling fields. The theory can then de- 
scribe a Yang- Yang anomaly [2] in which the second tem- 
perature derivative of the chemical potential along the 
phase boundary, namely, (fpa/dT"^, diverges at the crit- 
ical point with the specific heat exponent a. This is 
shown to entail, also, a leading term in the coexis- 
tence curve diameter that dominates over the previously 
known term [15]. The strength of the Yang- Yang 

anomaly, TZ^, is directly related to the linear pressure- 
mixing coefficient in the ordering field h that we have 
labeled j2 • 

For convenience of reference, we identify here the main 

definitions introduced and the explicit results derived. 
The coefficients ji, j2, ko, n^, Vi, V2 entering the 
nonlinear scaling fields to quadratic order are defined 
via (1.4)-(1.7), in terms of reduced pressure, p, chemi- 
cal potential, p and reduced temperature deviation t = 
(T — Tc)/Tc: see (1.7). The scaling ansatz, in the form 
developed is presented in (2.1)-(2.3) and the basic prop- 
erties of the overall scaling fmiction W±{y, j/4, 2/5, ■ • •) are 
set out in (2.4) together with (2.8) and (2.10) for small 
argument, \y\ 0, and in (2.12) for \y\ 00. Explicit 
expansions in powers of t for the phase boundaries, Pa{T) 
and pij{T), are reported in (3.11) and (3.15). The cor- 
responding expressions for the coexistence curve (and its 
diameter), and for the entropies of coexisting gas and 
liquid are (3.20) and (3.21); Eq. (3.31) gives the entropy 
on the critical isochore, p — pc- The linear mixing co- 
efficients, themselves are related to measurable critical 
amplitudes, etc. in (3.22) and (3.32)-(3.35). 

Section IV addresses the behavior of various critical 
loci in the p, p, and p vs. T planes. The critical isokyme, 
defined by ^ = ^c, is described in (4.5) and (4.7); the crit- 
ical isobar, in (4.13) and (4.15); the isotherm, in (4.19) 
and (4.21); and the critical isochore in (4.28) and (4.30). 
The family of k- susceptibility loci is defined via (4.31) 
while (4.38), (4.41) and (4.46), with (4.34) elucidate their 
behavior in the p, p, and density vs. T planes, respec- 
tively: see Figs. 1-3. The analogous /c-heat-capacity loci 
are defined in (4.50) while (4.62) describes their appear- 
ance in the (p, T) plane: sec Fig. 4. For both the k- 
susccptibility and /c-heat-capacity loci there are "opti- 
mal" values of k. depending only on the pressure- mixing 
coefficient j2, for which the dominant singularity at crit- 
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icality vanishes. Thus these particular loci approach the 
critical point in the (p, T) plane almost linearly: see 
(4.37) and (4.64). 

The magnitudes and inter-relations of the various crit- 
ical loci, including analytically continued loci, in the 
{p, T), (/i, T) and (p, T) planes are illustrated in Figs. 5, 6 
and 7, respectively. Corresponding analytic results for a 
general classical fluid, defined via the Landau expansion 
(5.4), are given in (5.6)-(5.9) and (5.10)-(5.19) and, for 
the A;-loci, in (5.22). The Yang- Yang relation and its gen- 
eralization to a linear density locus [in (5.34) and (5.40)] 
are discussed for classical systems for which there is 
no uniquely defined "anomaly" in Section V.D. Fig. 
8 illustrates that both pressure and chemical potential 
variation contribute to the specific-heat jump in a van 
der Waals fluid. 

Finally, in part II of this article [29] , the present scal- 
ing formulation with pressure mixing will be extended to 
finite-size systems (with periodic boundary conditions). 
The results, which include the definition and analysis 
of "Q-loci" , are of practical importance in the analysis 
of simulation data for near-critical asymmetric systems 
[6,25] as will be further demonstrated [29]. 

ACKNOWLEDGMENTS 

Erik Luijten collaborated in producing the coexistence 
curve estimates for the restricted primitive model elec- 
trolyte presented in Fig. 2. Support from the National 
Science Foundation (through Grant No. CHE 99-81772) 
is gratefully acknowledged. 



PcMs = ai3 + [a3i(a3i«4o - 4a2ia3ia|oa50 

- 40400,22 + 40,210,1)0,41 -|- 5031040050 

— 3a2i04oa6o) — 4032021040 + 050(4031043041 

— 4021022040 + 2O21O50 — 3O21O40O60) 

- 2051021040 +070021040] /8a|o- (A7) 
Pc = PcMo - ooo, Pi = PcMi - aoi, (A8) 

P2 = - O02 + 021/4040, (A9) 

P3 = - ao3 + 021(4022040 - 2O21O40O41 

-|- 021040 060 — O31O40 + 2O21O31O40O50 

- oii0^o)/8«40- (AlO) 

To specify the corresponding coefficients Oj^ for the 
van der Waals equation [see K(App. C)] we take Ac to 
be the thermal de Broglie wavelength at T = in a 
rf-dimensional system and define 



A = ln(AJf/26) - 1. (All) 
The leading nonvanishing coefficients are then 



o-oo 


= (|A-3)pc, ooi 


— doo 


-Pc, - 


-002 — 3ao3 - 


(A12) 


OlO 


= Ooo +Pc, On = 


fApc, 


-O12 = 


= 3oi3 = 2pc, 


(A13) 


a2o 


= 0, 021 = 3pc, 


040 = 


= O41 = 




(A14) 


^50 


3 

= ~^Pc, 0-60 = 


80 ^^C' 


O70 = 




(A15) 
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